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Abstract 

Current equation-of-state formulations, used for finite-difference cratering 
flow calculations, are cast into a framework permitting comparison of peak 
pressures attained upon impact of a sphere, with a half-space, along the Impact 
symmetry axis, to one-dimensional Impedance match solutions. On the basis of 
this formulation and application of thermochemical data, the regimes of melting 
and vaporization are examined. For the purpose of identifying material which 
will, upon isentroplc release from the impact- induced shock state, result in a 
solid just brought to its melting point, i.e., inclpiently melted (IM) ; com- 
pletely melted (CM); Just brought to its boiling point, i.e., Inclpiently va- 
porized (IV) ; and completely vaporized (CV) state, the pressures at which the 
critical isentropes intersect the Hugoniots of iron and gabbroic anorthosite 
(GA) are examined in detail. The latter rock type is assumed to be representa- 
tive of the lunar highlands. The Hugoniot pressures, for which IM, CM, IV, and 
CV will occur upon isentropic expansion, are calculated to be 2.2, 2.6, 4.2, and 
16,8 Mbar, respectively. For the high pressure phase (hpp) assemblage of GA, 
modelled as a mixture of plagioclase in the hollandite structure and pyroxene 
in the perovsklte structure, IM, CM, IV, and CV are calculated to occur upon 
isentropic expansion from Hugoniot states at 0.43, 0.52, 1.02, and 5.9 Mbars, 
respectively. The spatial attenuation of shock pressure along the impact axis 
is found to be clearly represented by two regimes, if the peak pressure, P, and 
radius normalized to that of the projectile, r, are fitted to expressions of the 
form Pt^r . At distances from 1.3 to 2.1 projectile radii into a GA target, the 
constant, a, is on the order of -0.2. This low-attenuation rate, near-field 
regime, extends further into the target at the slower Impact velocities and 
arises because of the slightly divergent flow associated with the penetration 
of a spherical projectile. For the near-field impact regime, an impact at 
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5 km/sec of an iron object v;ith a GA surface will induce CM for GA but the iron 
will remain solid. At 15 km/sec, partial vaporization (PV) occurs for both GA 
and iron, whereas at 45 kra/sec, CV occurs in both materials. In the case of a 
GA meteoroid striking a GA surface at 5 km/sec, partial melting results at 5 km/sec, 

PV occurs at 15 km/sec, and CV occurs in the near-field at 45 km/sec. At greater 
radii, in the far-field regime, the exponent, a, varies systematically from -1.45 
to -2.15 for Impacts of GA onto GA as the impact velocity is Increased from 5 to 
45 km/sec. For an iron projectile, the exponent, a, varies from -1.67 to -2.95 
as the impactor velocity increases from 5 to 45 km/sec. By comparison, the 
equivalent value of a, reported for both contained and surface explosions in 
various rocks is - 2 . It is suggested that, given field data on shock attenua- 
tion (based on identification of various shock metamorphic features versus distance) , 
overall crater size, and some chemical data as to the type of meteoroid which 
produced a crater, quantitative bounds on the Impact velocity of the meteorite 
may be obtained. 

Introduction 

The formulation of an equatlon-of-state for a rock type which we believe is 
typical of the lunar highland province, gabbroic anorthosite (GA) , provides an 
opportunity to investigate two important and related problems regarding Impact 
phenomena. These are: 

1) The impact velocities, and hence shock pressures, required to bring the 
meteorite and target material to the melting point, i.e.. Incipient melting (IM) , induce 
complete melting (CM), and produce a liquid at the boiling point, i.e., incipient 
vaporization (IV), and completely vaporize (CV) upon isentropic release is examined. 

We consider a highland terrane cratered by meteorites having equatlons-of-state 


ranging from iron to GA. 
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2) The spatial attenuation of peak shock-pressures In a GA composition, 
again, as the result of Impact of objects having equatlons-of-state similar to 
Iron and GA. 

It has been only recently determined that, even In very large and ancient 
astroblemes, geochemical evidence (via analyses of the minor siderophile element 
contents of Impact melts and ejecta ( Morgan et al. . 1975; Lambert , 1976)) allows 
specification of the causative meteorite. Moreover, recent quantitative studies 
of the effect of shock on minerals and their application to the spatial distribu- 
tion of shock metamorphic features within in-sltu shocked rock, can be used to 
Infer distinct bounds on the spatial shock attenuation rate ( Robertsen and Grieve , 
1977). If in addition the spatial attenuation rate can be related to the meteorite 
size and velocity, then improved estimates can be obtained of the chemistry, mass, 
and velocity of the meteoroids which have bombarded the earth's zone as a func- 
tion of time (Pence . 1972). Much of the required data Is not now available for 
lunar craters. However, if the moon formed and evolved near the earth, the late 
cratering history of the earth should, with minor modification, be applicable 
to the moon. 

In the present paper we obtain, within the framework of a two-phase Tlllotson 
equation-of-state model the impedance-match solutions specifying the peak pressure 
upon impact of a GA lunar surface by meteoroids of GA and iron in the 5 to 45 
km/sec range (O'Keefe and Ahrens , 1975, 1976). We then apply the available 
thermochemical data for the low pressure phase (Ipp) assemblage and a theoretical 
model for the high pressure phase (hpp) assemblage to construct critical release 
adlabats which define regimes in the pressure (P) - volume (V) plane which specify 
IM, CM, IV, and CV. Finally we fit the results of a series of hypervelocity 
impact, finite-difference, flow calculations specifying the spatial on-axis 
peak pressure attenuation due to the shock wave and compare these to previous 
calculations and measurements of pressure attenuation from explosions. 
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Conditions for One-Dlmenslonal Impact 

It is useful for the purposes of (a) providing an analytic comparison with 
finite-difference impact calculation and (b) obtaining an estimate 
of the impact velocity required to induce melting and vaporization in the target 
and impactor, to obtain one-dimensional Impedance match solutions for an assumed 
equation-of-state. As in our previous treatments of impact flows produced upon 
the impact of a meteorite with a planetary surface ( (O’Keefe and Ahrens , 1975, 1976), 
we employ a P-V-E (specific internal energy) , dependent equation-of-state of the 
form first proposed by Tlllotson (1962) . In the compressed region this is given 
by 

I + + BP^ (1) 


where a = 0.5 is the polytropic constant minus 1, at high temperature. The 

3P 

constant, b, is defined such (a + b) is the STP Griineisen parameter, y = V — , 

and A is the bulk modulus. Here r) = V /V and y = n-1, where V is the STP 

o o 

specific volume. The parameters b, B, and are obtained by fitting to Thomas- 
Fermi calculations of the equation-of-state in the 10^ 1-Ibar range, and Mugoniot data. 
The parameters for GA (Table 1) are based on previous shock wave measurements on 
lunar sample 15418 (Ahrens et al. , 1973). As in O’Keefe and Ahrens (1976), we 
consider the equilibrium equation-of-state for GA to have discreetly different 
Tlllotson parameters which specify the Ipp and hpp assemblage. For the latter, 
hpp, assemblage, stable above ''>150 kbar, we assume plagioclase has the density 
corresponding to the hollandite structure and pyroxene has a density corresponding 
to the perovskite structure. The effect of applying this relatively complex, 

REPRODUCIBILITY OP THE 
QEIGINAL PAGE IS POOR 
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but realistic equation-of-state to the description of the impacts onto GA is 

examined in some detail by O'Keefe and Ahrens (1976). 

For convenience we define, E = 0, for the Ipp at STP. To invert Equation 1 

to the P“U (particle velocity) plane, the energy term for a Hugoniot state in 
s 

the hpp regime is given by 

E = - V)/2 - Ejj (2) 

where is the Ipp, specific volume and E^^^ is the Increase in internal energy 
in going from the Ipp to the hpp assemblage at STP. 

It should be noted that if the E^^^ term is not included in Equation 2, the 
resultant Hugoniot curve generated would be that of the metastable Hugoniot 
( McQueen et al. , 1967) . Moreover we note that in applying the definition of n 
and M for describing the hpp, the initial specific volume used, V^, is that of 
the hpp. In calculating the Hugoniot curve for GA (Ipp) and for iron, 
and Upon eliminating E between Equations 1 and 2, the Hugoniot pressure, 

P = Pjj, is given analytically upon solution of the following quadratic equation 
where V (and hence p and also n) is the independent variable 

a'P^ + b'P„ + c’ = 0 (3) 

H H 

Here 

, s I as~l 1 

^ ■ k [ V J 

b' = (a + b)s/V -1+ E+R(l-2as/V)/k + As/k 
c' - a/(kV) - Ej^(a + b)/V + 2(l-Ej,^/k) 

where 


(4a) 

(4b) 

(4c) 
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s » (V - V)/2 
oo 

k = E n ^ 

o 

£ == IjA + 


(4e) 

(4f) 

(4g) 


The negative sign is used in the usual quadratic formula. Once P|j(V) along the 
Hugoniot, centered at is calculated, the corresponding value of u^ is simply 

calculated from the Rankine-Hugoniot equation 


By solving Equation 3 at a series of specific volumes, the relations given 

in Figure 1 are obtained from the parameters in Table 1. Straightforward appli- 
cation of the method of Impedance matching (Duvall and Fowles , 1963) (matching 
press'ire and particle velocity), and application of the Rankine-Hugoniot equa- 
tions yield the peak shock states defined in Tables 2 and 3 for Impact of iron 
and GA projectiles on a GA half-space. 


To calculate the amount of melting and vaporization that occurs in an impact 
event, the release isentropes for IM, CM, IV, and CV were calculated using the 
Tillotson equation-of-state (Equation 1). The calculation of the isentropes 
requires that at least one state along the desired isentrope be known initially. 
These initial states are typically the internal energy or temperature at one 


atmosphere pressure. To obtain the internal energy contents, at one atmosphere 
pressure of iron required for IM (1809°K) , CM, IV, (3145°K) , and CV, we have 
Utilized enthalpy and entropy data tabulated in the JANAF Tables (1965). In the 
case of GA, the above melting and vaporization processes are Incongruent, and an 
approximate treatment was used for simplicity. Following the rock modelling scheme 
of Ahrens and O'Keefe (1972), GA was assumed to comprise a mixture of 0.714 (mass 
fraction) anorthite, and 0.286, enstatlte. Using the JANAF thermo chemical data. 
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through the MgSi 03 melting point (1798“K) and Into the liquid, and extrapolating 
using the entropy systematlcs described in Ahrens and O'Keefe (1972), the energy 
and entropy for incipient and complete vaporization is inferred. Assumptions 
regarding volitllzatlon are made on the basis of oxide volatility and are dis~ 
cussed by Ahrens and O'Keefe (1972). In the case of CaAl2Sl208, a similar but 
even less certain extrapolation was employed using the thermo chemical data tabu- 
lated in Robie and Waldbaum (1968) and, again, the systematlcs of Ahrens and 
O'Keefe. The densities, at one atmosphere, given in Table 1 for IM, CM, and IV for 
the Ipp and hpp assemblages of GA and iron are obtained from Equation 1. 

Only in the case of iron can meaningful comparisons be made with other 
theoretical treatments and the actual densities at CM and IV. Table 1 gives 
densities of 7.36 and 6.10 g/cm^ for CM and IV, respectively. Interpolating, 
the results of the Mle-Griineisen reduction of the iron data by McQueen et al. 

(1970) to 1809“ yields a considerably lower density, 6.696 g/cm^, for CM. The 
CRC Handbook (1970) lists density values ranging from 6.99 to 7.24 g/cm^ for 
iron at the llquidus, slightly lower, but in good agreement with those calculated 
from Equation 1. For IV, at 3145“K, the CRC Handbook gives density values ranging 
from 5.84 to 5.93 g/cm^ which again is close to, but still slightly lower than, 
the 6.10 g/cm^ value indicated in Table 1. 

The onset of IM, CM, IV, and CV upon expansion (assumed isentroplc) from 
Hugoniot states even for the one-dimensional case is of interest. We have 
calculated the P-V isentropes which pass through the states IM, CM, and IV 
specified in Table 1 at atmospheric pressure and high temperatures using Equation 
1 and found their intersection with the Hugoniot curves (Table 4). In the case 
of iron and the Ipp of GA, this is a straightforward procedure which merely 
involves numerically calculating the locus of V, P, and (energy) states reached 
by an isentropic process which satisfies the energy integral along the isentrope 
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Es = - I V(P,E) dP 
•'V(IM,CM,orIV) 

In principle, the actual value of the entropy along one of the Isentropes, so 
calculated, is not explicitly specified but, in fact, requires an additional 
assumption equivalent to a thermodynamic model for CpOr C^, the specific heat at 
constant pressure or volume. Instead, for the Ipp we have used the available 
thermodynamic data for enstatite and anorthite, to separately estimate the 
entropy associated with the critical isentropes for IM, CM, and IV using the 
method outlined in Ahrens and O'Keefe (1972). (These data are utilized In calcu- 
lating the critical Isentropes for the hpp assemblage below.) 

In the case of CV, the pressure is obtained by using the distended formula- 
tion of the Tillotson equatlon-of-state: 



! bE/V 

1 - 2+1 


^^-e[(v/v^-i)-i] 




( 7 ) 


in* the regime where 

1 > V/V for E > E^„ 

o CV 

where a = 3 = 5, and is the CV energy at standard pressure. In the partial 
vaporization regime 


Vo > V and < K < 

where E^^ is the energy associated with IV. We utilize the interpolation relation 
( Allen , 1967; Hageman and Walsh . 1970) 
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(E > E^.^) P„ + (E^„ - E) P 


CV 




( 8 ) 


where P_ refers to the pressure calculated from Equation 7, the distended region, 
r« 

and Pp refers to the pressure calculated from Equation 1, the compressed region. 

The straightforward calculation of the Isentrope (Equation 6) using the 
distended formulation Tlllotson equatlon-of-state and starting at a one atmosphere 
pressure state at complete vaporization leads to unreasonable pressures (2^ 10^ Mbar) 
for the Intersection of the complete vaporization isentrope with the Hugoniot. 

The reason for this is that the second term In Equation 7 does not have the 
correct form in the distended regime, where u Is negative. The correct cohesive 


energy, _ Is not obtained by taking 


/ ysco 

V=Viv 


. This problem was partially 


circumvented by choosing a reference state for complete vaporization, 
which did not involve the distended region of the equation-of-state. 

To this end, we adopted the approach of Zeldovich and Ralzer (1967), who demon- 
strate, for single phase materials, that the isentrope that passes through the vapor 
liquid critical point has a specific internal energy equal to approximately twice 
the binding energy at standard volume. Because of this assumption the curves 
shown In Figures 2 through 4 for CV should be considered approximate. 

The construction of release isentropes for the hpp assemblage (Figure 4) 
which is used as a model for GA, requires two independent assumptions in addition 
to the parameters of the P, V, E equation-of-state: 

1) The entropy difference between the Ipp and hpp assemblage; no data 
presently exists giving the slope of pertinent phase boundaries. 

2) A model for the specific heat for the assemblage. 

Although more elaborate methods have recently been suggested for estimating 
entropies of phases based on detailed knowledge of crystal structure (e.g. , Saxena . 
1976), we have used the method of Fyfe (1958) to obtain the entropy of the hpp 
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assamblage. This method sums the entropies of the equivalent oxide mixture 
(data obtained from Roble and Waldbaum , 1968). By using the STP entropy of 
stlshovlte, this procedure probably provides a good approximation, as the major 
effect of the assumed phase changes in both pyroxene and plagloclase is the 
Increase of Si"^** coordination from four to six 0“^ Ions. Upon mass weighting 
entropies for the plagloclase and pyroxene hpp mixture, the transition entropy 
at STP is given by; 

oqooif 

Ipp hpp ; = -0.188 X 1Q7 ergs/gm“K (9) 


Recognizing the shortcomings of the Debye theory for describing C^, when 
applied to silicates (Kieffer and Kamb , in preparation) , we have assumed an 
equivalent STP, Debye temperature, 0^, for the hpp assemblage. Using the density- 
Debye temperature relation proposed for oxides by Anderson (1965) , a value of 
0JJ = 1029®K is obtained. The entropy, at high temperature at atmospheric pressure 
for the hpp assemblage. Is then calculated from 



where, taking into account the entropy of the known Ipp assemblage, and Equation 

9, = 0.534 X 10^ ergs/gm°K. 

hpp 

From the thermodynamic relations: 



/[“'■&) (lib’] 


( 11 ) 
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and from Equation 1, it follows that 


/av\ ( 

Wp ■ \ 



bE 


(12) 


and 


. b bE 
Y “ a + 2 

X X 




(13) 


where x “ E/E 


The variation of 0^ with volume follows from 


(90p/8T)p = - ya 


(14) 


where 


a = 


(s) [‘'Ks) "] 


(15) 


In order to obtain the internal energies, relative to the Ipp at STP, for the 
hpp at one atmosphere, such that the entropy value associated with the critical 
isentropes for IM, CM, IV, and CV are the same as for the Ipp, we have simul- 
taneously solved Equations 10-15, using T, the absolute temperature as the in- 
dependent variable. Thus upon isentropic release from a shock state either in 
the Ipp or hpp regime to one atmosphere and high temperature, and achieving 
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states corresponding to IM, CM, IV, or CV, the entropy of the material is the 
same, however, the Isentropi paths from high pressure differ for the Ipp and 
hpp cases, inferring that reverse transformation from the hpp occurs. 

We have not yet compared the results of the Debye theory, when used with 
the Tillotson formulation as described in Equations 12-15, when the high 
temperature entropy has been explicity measured. 

To describe the low cohesion of iron or GA, we have assumed P^, = 0, when 
V^/V < 0.99 and V^/V < 0.995 respectively. This corresponds to a dynamic tensile 
strength of 0.64 and 0.71 kbar, respectively. The latter physical quantities 
are both poorly constrained experimentally, although recent measurements of the 
tensile strength of novaculite and quartzites ( Shockley et al. . 1973) report 
some data in this range. By assuming = 0 upon tensile failure, a physical 
model of the formation of a fragment dust and coexisting vapor is implicitly 
assumed. The total pressure upon failure results from only the gaseous (P^.) 
component of a mixture of dust and gas. Using Equations 1, 6, and 7, and the 
assumption concerning dynamic brittle failure to dust, the release isentropes 
shown in Figures 2, 3, and 4 are calculated. 

Results 

Applying the above detailed formulation of the equation-of-state of GA, the 
impact flow nd energy partitioning resulting from the interaction of a series 
of hypothetical iron and GA meteorites with a GA planetary surface (half-space) 
are calculated for impact velocities ranging from 5 to 45 km/sec. Previous 
results for energy partitioning and ejecta distributions upon impact of a 10 cm 
diameter, 15 km/sec iron object are given in O’Keefe and Ahrens (1976). 

Ejecta (0.5 cm, initial zone size) distributions and energy partitioning 
versus impact velocity will be presented elsewhere (O'Keefe and Ahrens, 
in preparation). Because the shape of the isentropes for the hpp 
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assemblage of GA strongly controls peak pressure decay, c.f., O'Keefe 
and Ahrens (1976), and field data relating to spatial shock attenuation are 
beginning to become available, c.f., this volume ( Robertsen and Grieve , 1977), 
initial calculations of shock attenuation along the impact, symmetry axis are 
presented in Figures 5 and 6 and Table 5. The scatter is largely the result 
of the rezoning procedure in the calculation. As previous workers have noted 
(e.g., Gault and Heitowit , 1963; Bjork and Rosenblatt , 1965; Heyda and Riney , 1965), 
two regimes of attenuation are to be expected. These are indicated as near- and 
far-fleld in Table 5. Within the precision of the present calculations, centerline 
peak pressure decay in the near- and far-field is described in terms of: 

logfoCP - Mbar) - a log^^CR/R^) + b (15) 

where P is the centerline pressure, R is the distance from the point of impact and 
the initial meteoroid radius. In the near-field regime, the peak pressure 
may be obtained from the one-dimensional impedance match solution (Tables 3 and 4) . 
The agreement with impedance match solutions indicated in Figures 5 and 6 appears 
to be closer, the higher the impact speed, as R ->• 0. This may be due to shorter 
time steps in the calculation at the higher impact speeds. The slow spatial decay 
rate observed in the near-field is insensitive to both impact velocity and 
impactor type. In the present case of a spherical projectile, the near-field 
shock attenuation undoubtedly results from the slightly divergent flow induced 
in the target as the projectile imbeds itself into the target. In the case of 
the normal impact of a flat-nosed cylinder or rod, or in the assumed constant- 
hemisphere energy model of Gault and Heitowit, no attenuation, i.e., a = 0, is 
expected in the near-fleld regime. The onset of the characteristic far-field 
attenuation rate begins at depths corresponding to >1.6 to '\-2.1 projectile radii 
for the case of , iron impacting GA as the impact velocity decreases from 45 to 
5 km/sec. In the case of GA impacting GA, the far-fleld rates are established 
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at depths of from ''>1.3 to 'V/1.6 projectile radii, as the Impact velocity decreases 
from 45 to 5 km/sec. Dienes and Walsh (1970) have shown as part of their late- 
stage equivalence studies of hypervelocity impact that late-stage far-field flows 
are insensitive to the exact shape of the impactor, provided the axial to radial 
dimensions of the impactor are roughly comparable. 

Comparison of the present results v/ith other attenuation calculations and 
measurements are of interest. Dienes and Walsh found good agreement between 
their computer solution with experimental data, over a range of 'vlQ to 1000 kbar, 
for the impact at 7.3 kra/sec of an aluminum sphere with an aluminum half-space. 

The centerline particle velocity in this case, which is to first-order propor- 
tional to shock pressure, decays as, a = -2, which is comparable to che value 
of "a", listed in Table 5, for the far-field, in this impact velocity range. A 
larger set of experimental impact data, up to km/ sec, for various aluminum 
projectiles is shown to satisfy a relation comparable to present results, where 
a = -1.6 (Billingsley , 1969). Gault and Heitowit (1963), on the basis of their 
constant energy, hemisphere model, calculate, for the impact at 6.35 km/sec of 
an aluminum sphere onto basalt, (similar to the present An ->• An calculation) an 
initial pressure decay rate with a = -3 to -4, which is considerably greater 
than predicted here. However, at the point where the shock has decayed to pres- 
sures on the order of '^>10^ kbar, their indicated attenuation rate is considerably 
slower than is predicted here. 

Since the cratering effects of large explosions are often used for scaling 
of impact flows to greater dimensions than those available from results of experi- 
ments, it is of interest to compare pressure decay rates from near-surface 
and contained explosions. A large number of such data are summarized by 
Cooper (1973), who finds a decay coefficient of the centerline pressure of a 
= -2 is compatible with field data for both contained and surface explosions. 
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Thl3 is close to the value found, for example, for An -»• An at 15 km/sec. Also, 
Butkovlch and Borg (1974) demonstrated that the decay of the experimentally 
observed shock pressure In the range 650 to '^l kbar surrounding the 5 kton, 

Hardhat nuclear explosion in granodiorite also gave a a = -2 dependence. More- 
over, they point out that in general, calculatlonal results such as presented 
here, always give, to a first approximation, the decay rates in terms of pre-flow, 
l.e., Eulerian coordi^ Ttes, whereas attenuation inferred from post-flow shock 
metamorphlc features, must for meaningful comparisons, be corrected such that 
material motion, after the shock wave has passed, is taken into account. 

Discussion 

The major new result described here, is that different impact velocities, 
and meteorite shock impedances give rise to significantly different spatial peak 
shock attenuation rates in the far-field. 

We anticipate that present calculatlonal results will prove useful in placing 
bounds on the shock pressures, experienced by both In-situ rocks and ejecta, in 
the vicinity of impact craters, which petrologic and geochemical analyses indicate 
are partially melted, melted, or, Incongruently vaporized. Moreover, (a) given the 
size of a crater, and hence a measure of the total projectile energy, (b) some 
knowledge of the chemistry, and hence shock impedance, of the impactor and (c) a 
measure of the spatial rate of peak shock-pressure attenuation (probably Inferred 
from shock metamorphlc features) we believe it will be possible, on the basis 
of the present results to estimate the Impact velocity associated with a given 
crater. Although the required knowledge of crater size, meteoroid impedance, 
and spatial shock decay rate is now limited to only a few terrestrial craters, 
we hope the methodology presented will be useful in inferring impact 
velocities and cratering histories for other terrestrial planets. 
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Table 2 


Peak Shock (Impedance Match) States for Gabbrolc Anorthosite 
Impacting Gabbrolc Anorthosite 


Impact 

Velocity 

(km/sec) 

Shock 

Pressure 

(M)ar) 

Particle 

Velocity 

(km/sec) 

Shock 

Velocity 

(km/sec) 

Density 

(gm/cm^) 

Internal a) 
Energy Density 

(10^® ergs/gm) 

5 

0.62^ 

2.50 

8.45 

4.17 

3.13 

7.5 

0.99^ 

3.75 

8.99 

5.04 

7.03 

15.0 

3.04*^ 

7.50 

13.81 

6.43 

28.13 

30.0 

10.11*^ 

15.0 

22.96 

8.47 

112.50 

45.0 

21.29^ 

22.50 

32.23 

9.73 

253.13 


a) With respect to low-pressure phase of gabbrolc anorthosite 
at STP. 

b) Partial melting upon isentropi.c release. 

c) Partial vaporization upon isentropic release. 

d) Con?)lete vaporization upon isentropic release. 


Peak Shock (In?>edance Match) States for Iron Inqpactiixig Gabbrolc Anorthosite 
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Table 4 


HugonioC Pressures (Mbar) Inducing a Change of 
State Upon Isent topic Release to Atmospheric Pressure. 


Iron 

Gabbroic 

Anorthosite 


(^PP) 

(hpp) 


Incipient 

Melting 

2.2 

0.97 

0.43 

Complete 

Melting 

2.6 

1.1 

0.52 

Incipient 

Vaporization 

4.2 

1.9 

1.02 

Complete 

Vaporization 

16.8 

7.9 

5.9 



Table 5 


Peak Centerline Pressure Attenuation; 
logl 0 (P“Mbar) = a logig(R/R^) + b 


Near-Field Far-Field 


An An 

a 

b 

r2<^> 

a 

b 

2 

5 km/sec 


— 

— 

-1.45 

0.151 

0.85 

15 km/sec 

-0.222 

0.285 

0.83 

-1.97 

1.399 

0.95 

45 km/sec 

-0.293 

1.208 

0.92 

-2.15 

2.373 

0.92 

Fe An 

5 km/sec 

-0.232 

-0.398 

0.69 

-1.67 

0.681 

0.91 

15 km/sec 

-0.295 

0.470 

0.88 

-2.49 

1.475 

0.97 

45 km/sec 

-0.149 

1.220 

0.89 

-2.95 

2.082 

0.99 


(a) Linear regression correlation coefficient. 



'igure Captions 


1. Pressure-particle velocity Hugonlot relation for iron, and low-pressure 
phase (£pp) and high-pressure phase (hpp) assemblage of gabbrolc 
anorthosite. Gabbrolc anorthosite Hugoniots are both centered at spe- 
cific volume of Jlpp and calculated using the Tlllotson parameters of 
Table 1. 

2. Hugonlot (H) for iron and calculated critical Isentropes, corresponding 
to incipient melting (IM), complete melting (CM), incipient vaporization 
(IV) and complete vaporization (CV). The H and IM curves agree 
closely with similar results obtained by Tlllotson (1962). 

3. Hugonlot and critical Isentropes for £pp, gabbrolc anorthosite. Curve 
labeling similar to that in figure 2. 

4. Hugonlot and critical Isentropes for hpp, gabbrolc anorthosite. Curve 
labeling similar to that in Figure 2. Conq>ression, V/V^, corresponds to 
zero-pressure specific volume of hpp. 

5. Logj^Q peak shock pressure, versus, normalized radius at various inq>act 

velocities. is radius of in5>actor for iron object. Parameters of 

lines fit through calculation are given in Table 5. Arrows (IMS) 
indicate impedance-match solutions (Table 3) valid for one-dimensional 
flow. 

6. peak shock pressure, versus, normalized radius at various Impact 
velocities for gabbrolc anorthosite impactors. Arrows Indicate one- 
dimensional flow pressures (Table 2), Notation is similar to Figure 5. 
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